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Abstract 

Some properties of cosmological models with a time variable bulk 
viscous coefficient in the presence of adiabatic matter creation and 
G, c, A variables are investigated in the framework of flat FRW line 
element. We trivially find a set of solutions through Dimensional 
Analysis. In all the studied cases it is found that the the behaviour of 
these "constants" is inversely proportional to the cosmic time. 



1 Introduction. 

Recently several models with FRW metric, where "constants" G and Lambda 
are considered as dependent functions of time t have been studied. For these 
models, whose energy-momentum tensor describes a perfect fluid, it was 
demonstrated that G oc t a , where a represents a certain positive constant 
that depends on the state equation imposed while A oc t~ 2 is indepedent of 
the state equation (see Q). More recently this type of model was generalized 
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by Arbab (see P]) who considers a fluid with bulk viscosity (or second viscos- 
ity in the nomenclature of Landau (see 0)). The role played by the viscosity 
and the consequent dissipative mechanism in cosmology has been studied by 
many authors. The heat represented by the great entropy per baryon in the 
cosmic background radiation provides a good indication of the early universe 
and a possible explanation of the huge entropy per baryon which is believed 
to have been generated by physical dissipative processes acting at the be- 
ginning of evolution. These dissipative processes may be responsible for the 
smoothing out of initial anisotropies. 

In the models studied by Arbab constants G and A are substituted by 
scalar functions that depend on time t. The state equation that governs the 
bulk viscosity is: £ oc £oP 7 where 7 is a certain indeterminate constant for 
the time being 7 e [0, 1]. Amongst all the possible 7 we are only going to 
concentrate on 7= 1/2 since, as pointed out by Golda et al. (see M) upon 
demonstrating that for an adequate election of the state equation the viscous 
models are topologically equivalent (structurally stable) to the classic FRW. 
Golda et al show that the bulk viscous fluids with FRW symmetries can 
structurally approximate the dynamics of the classic FRW. They found that 
approximation only takes place when the parameter of the bulk viscosity 
follows the state equation £ oc cjop 1 ^ 2 - They seek spatially homogeneous 
and isotropic solutions with bulk viscosity that could be approximated in 
a structurally stable way to the dynamics that describe the classic FRW 
described by an ordinary perfect fluid. More specifically, they looked for such 
spatially homogeneous and isotropic solutions which after being disturbed by 
the dissipative parameter have a dynamics that is topologically equivalent to 
the classic FRW models. Therefore the only viscous models that can be 
approximated to the classic FRW in a structurally stable way are those that 
follow a state equation £ oc £oP 1//2 because of their viscous parameter. 

The first authors that have studied the convenience of also considering the 
variation of " constant" c, the speed of light, as a scalar function depending on 
time t, to solve some of the problems that classic FRW models present, have 
been (see Others (see |7J and ||) demonstrate through a variational 

formulation that the resulting field equations continue being the Friedmann 
ones if a FRW metric is considered. In these papers it is not so much the 
accurate calculation of such a variation for c that is taken into consideration 
as the examination of the benefits of a modification in the field equations. 
Thus, these models are able to explain the problem of the horizon and the 
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flatness of the universe, in the same way as the inflationary models. 

Following Arbab's line, the author of this article (see 0) attempts to 
generalize the models with G, c and A variables through the consideration of 
a flow with bulk viscosity. In this type of models, (the structurally stable 
7 = 1/2) inspite of considering viscous fluids, that involve a natural pro- 
duction of entropy (irreversible processes), it demonstrates that the entropy 
continued being constant (as in the classic FRW). In this paper we consider 
both mechanisms of matter production and mechanisms of entropy. To at- 
tempt to solve this problem we consider in this paper so much mechanisms of 
matter production as of entropy. The problem of matter creation in the uni- 
verse has been studied by many authors following Prigogine and co-workers 
pioneering ideas (see JTUJ). We will follow Lima et al's work, (see ( [JTTJ] ) . This 
problem in the bulk viscous fluids framework has been studied by Desikan 
( [|TI| ) . The general idea is the following: to consider the universe as an open 
system, the creation of matter itself generatersemtropy and thus the second 
law of the thermodynamics can be incorporated into the field of equations. 
In the case of open systems the number of particles in a given volume is 
not fixed and though we consider that the transformation is adiabatic, it is 
demonstrated that the entropy is growing when the matter creation acts as 
an internal energy mechanism. In this paper we calculate the variation of 
" constants" G, c and A in the framework of a model with FRW symmetries, 
with k = i.e. flat, with adiabatic matter creation and bulk viscosity as sep- 
arate irreversible processes. We consider the processes of adiabatic matter 



creation (see [11]) to attempt to solve the problem of entropy. We wish to 
draw attention to how the use of the Dimensional Analysis (D.A.) (see ||13|| ) 
enables us to find the solution to our model in a trivial way. We do not want 
to carry out a meticulous study of the cosmological implications of the results 
obtained for each of the parameters on which the solution depends (lu, (3, 7) 
(postponing this paragraph for a subsequent article) but simply to demon- 
strate them. We will show that for an adequate election of such parameters 
we obtain the same behaviour for the principal quantities as those obtained 
for the classic FRW models. The behaviour obtained for the "constants" is 
that for these parameters they always vary in a way inversely proportional 
to the time. These results are compared to those of other authors. 

The paper is organised as follows: In the second section the governing 
equations of our model will be shown and considerations on the followed 
dimensional method will be made. In the third section we will make use of 
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the D.A. (Pi theorem) to obtain a solution to the principal magnitude that 
appears in the model and finally in the fourth section we will present two 
particular cases of the obtained solutions together with some conclusions. 



2 The model. 

For a flat universe k = with FRW symmetries i.e. we assumed homogeneity 
and isotropy and therefore there will not be spatial variations of the "con- 
stants" G, c and A solely temporary. We suppose equally that our fluid is 
bulk viscous (second viscosity) and we consider mechanisms of creation of 
matter. With these suppositions the equations that govern the model are as 
follows: 

J" if') 2 87cG(t). , jMi/ , 
" 2 y - KJ jf = -^±L(p + Pe ) + c\t)A(t) (1) 

n + 3nH - if> = (3) 

where n measures the particles number density, ip is the function that mea- 
sures the matter creation, H = f'/f represents the Hubble parameter (/ is 
the scale factor that appears in the metrics), p is the thermostatic pressure, 
p is energy density and p c is the pressure that generates the matter creation. 

The creation pressure p c depends on the function ip. For adiabatic matter 
creation this pressure takes the following form ( ||11| ): 

-p + p 



Pc 



-0 



(4) 



3nH 

The state equation that we will use is the known expression 

p = ujp (5) 

where u = const, to £ [0, 1] physically realistic equations, making in this way 
that the energy-momentum tensor verify the energy conditions. 

We need to know the exact form of the function ip the one which is de- 
termined from a more fundamental theory than involves quantum processes. 
We assumed that this function continues the following law: 

ip = WnH (6) 
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here we are following to Lima et al (see [II]) for other treatment |]12| while 
Prigogine et al [ 10 1 follows this other law ip = kH 2 ) where f3 is a dimensionless 



constant (if /3 = then there is matter creation since ip = 0) that it is given 
presumably by models of particles physics of matter creation. Physically 
one hope that the most interesting situations emerge during phases in those 
which i.e. (3 ~ 1 will be of the order of unity. 

The conservation principle carries us to the following expression: 

p' + 3(u, + l)p^ = (u, + l)/ (7) 
/ re- 
integrating the equation (|7]) we obtain the following relationship between 
energy density and the radius of the universe and what is more important 
the constant of integration necessary for our subsequent calculations: 

P = ^, /3 /- 3(w+1)(1 - /3) (8) 

where A^^ is the constant of integration that depends on the state equation 
that is considered i.e. of the constant u and of the constant j3 that measures 
the matter creation. 

The effect of the bulk viscosity in the equations is shown replacing p by 
p — 3£H where £ follow the law £ = £oP 7 (see fL4jl,|P] and f!2|l ). This last 



state equation, in our opinion, it does not verify the homogeneity principle 
by this reason we modify it by: 

£ = W (9) 

where the constant fc 7 causes that this equation yes will be dimensionally 
homogeneous for any value of 7. 

The dimensional analysis that we followed needs to make the following 
distinctions, we need to know beforehand the set of fundamental quantities 
together with that of unavoidable constants (in the nomenclature of Baren- 
blatt are designated as governing parameters). In this case the single one 
fundamental quantity that appears in the model is the cosmic time t as can 
be deduced with facility of the homogeneity and isotropy supposed for the 
model. The unavoidable constants of the model are the constant of inte- 
gration that depends on the state equation u and of the mechanisms 
of matter creation f3 and the constant fc 7 that controls the influence of the 
viscosity in the model. 
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In a previous work (||15||) was calculated the dimensional base of this type 
of models, being this B = {L, M, T, 9} where 9 represents the dimension 
of the temperature. The dimensional equation of each one of the governing 
parameters are: 

[t] = T [Awfi] = L 3 (-+ 1 )( 1 -«- 1 MT~ 2 

All the derived quantities or governed parameters in the nomenclature of 
Batrenblatt we will calculate them in function of these quantities (the gov- 
erning parameters) , that is to say, in function of the cosmic time t and of the 
two unavoidable constants fc 7 and Aj,/3 with respect to the dimensional base 
B = {L, M,T,9}. 

3 Solutions through D.A. 

We go to calculate through dimensional analysis D.A. i.e. applying the Pi 
Theorem, the variation of G(t) in function of t, the speed of light c(t), the 
energy density p(t), the radius of the universe f(t), the temperature 9 if), the 
entropy Sif) and the entropy density s(t),the viscosity coefficient £(t), the 
particle number density n(t) oc f~ 3 and finally the variation of the cosmo- 
logical "constant" A(£). 

The dimensional method carries us to (see 

3.1 Calculation of G{t) 

As we have indicated above, we are going to accomplish the calculation of the 
variation of G applying the Pi theorem. The quantities that we consider are: 
G = G(t, k y , Au t p). with respect to the dimensional base B = {L, M,T,9} . 
We know that [G] = L 3 M~ 1 T~' 2 

( G t k n \ 

L 3 7- 1 3(o; + l)(l-/3)-l 

M -1 1-7 1 

V T -2 1 2 7 - 1 -2 / 
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we obtain a single monomial that leads to the following expression for G 

2 2+3(^+l)(l-/3) 
r a 3(^ + l)(l-/3) 7 3(w + l)(l-/9)(-y-l) 

. r 2+3(^+l)(l-/3) I 
t- 4 -[3^ + l)(l-/3,( 7 -l)J 

3.2 Calculation of c(t) 

c(t) = c(t, kj, Aufi) where [c] = LT^ 1 =>- 

i i 

n (f\ 4 3(^ + l)(l-/3) , 3(^ + l)(l-/3)( 7 -l) 



1 [3(0, 



+ l)(l-/3)( 7 -l)J 



3.3 Calculation of energy density p(t) 

p = p(t, k 7 , Au^p) with respect to the dimensional base B, where [p] = 
L -i MT -2 

p oc kt^t^ (12) 

we observe that this relationship shows us that energy density does not de- 
pend neither on the state equation 00 nor on the mechanisms on creation of 
matter i.e. does not depend on the constant A^^ solely on the viscosity of 
the fluid. 



3.4 Calculation of the radius of the universe f(t). 

f = f(t, fc 7 , A Wi0 ) where [f] = L => 

f oc fc ff"+D(i- w-r-D f 3K1)( r- 1 ffl7 - 1) ( 13 ) 

We can observe that: 

- [?>(3{u + 1)( 7 - 1) - 3n(tv + 1) + 3u + 4] 

H _f'_ ( 1 \ 1 

/ " U W + l)(l-^)( 7 -l)J t 
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3.5 Calculation of the temperature 9{t). 

9 = 9(t, k y , A^^ks) where fee is the Bolztmann constant : [9] = 9 and 
[k B 0] = L 2 MT Li 

k B 9 oc A^ 1K1 ' m k ^™^ t~ t ( M+ D(i-^(-y-i) ] (14) 

3.6 Calculation of the entropy S(t) : 

S = s(t, fey, Au t /3, a) where a is the radiation constant. [S] = L 2 MT~ 2 9~ 1 



l-£( w +l)(l-0) 



S* oc A ( " +1)(1_/3) fcn" +1)(1_/3)(7_1) t [( w + 1 )( 1 -«(^- 1 )J a i (15) 

3.7 Calculation of the entropy density s(t) : 

s = s(t, fey, A^^, a) where a is the radiation constant, [s] = L~ 1 MT~ 2 9~ 1 



s oc A^kn^t'i^^ia* (16) 

3.8 Calculation of the viscosity coefficient : 

£= £(t, fey, 4^) where [f ] = L^MT" 1 

f oc feT^^r (17) 

3.9 Calculation of the cosmological constant: A(t). 

A = A(t, /c 7 , Aufi) where [A] = L~ 2 



-2 -2 2 



A oc A^5 +1H1 " W ^ ( " +1K1 ^ )( ^ 1) t 3(^ + D(i^)(,-i) (is) 

4 Different cases. 

All the following cases can be calculated without difficulty. But as we have 
indicate in the first section we are going to centre our attention only in those 
models that follow the law £ = k^p 1 / 2 i.e. 7 = (1/2) that corresponds to 
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models that are topologically equivalent to the classic FRW. We are going 
to study two models with 7= (1/2), one with uj = 1/3 that corresponds to 
a universe with radiation predominance and other with u = corresponding 
to a universe with matter predominance. 



4.1 7= 1/2 and u = 1/3 

7oc A^ 1 



1 1 



if /3 = G oc r 1 , c oc r 1 / 2 , A oc r 1 

c oc i -1 / 2 also has been obtained by Barrow and Troiskii ||, with respect 
to the rest of the quantities we have obtain the same behaviour that Lima et 



al. |TT 



p oc k^t 2 p oc t 

4(1-/3) 7, 2 ~2(l-^i) ,-2H 



-2 



f OcA^kn^-'H^ 

1 3 3/3 3fl 

aTSoc^ (1 -^„ 2M t5(H3) 

_! 3 3 _3 

a~ S OC A^knt~2 soct 

e oc ^r 1 e cx t- 1 

Si /? = 6 oc r 1 ' 2 , / oc t 1/2 , 5 oc t° = const. 

With respect to the thermodynamic behavior, the matter creation formula- 
tion considered here is a clear consequence of the nonequilibrium thermody- 
namic in presence of a gravitational field. We see that the j3 parameter works 
in the opposite sense to the expansion, that is, reducing the cooling rate with 
respect to the case where there is no matter creation. A very meaningful re- 
sult is the fact that the spectrum of this radiation cannot be distinguished 
from the usual blackbody spectrum at the present epoch (see ||11|| ). Therefore 
models with adiabatic matter creation can be compatible with the isotropy 
currently observed in the spectral distribution of the background radiation. 
We observe equally, that the obtained model is clearly irreversible (classic 
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FRW is reversible). We want also to express the fact that all the important 
thermodynamic quantities of the classic FRW models are recovered if we 
made (3 = 0. (see JTB|) 



/ oc t 1/2 , p oc t~ 2 , 9 oc r 1/2 , S = const., s oc r 3/2 

Finally it is interesting to stick out that the presented model may signifi- 
cantly alter the predictions that make the classic FRW on the abundance of 
elements. Such result puts a possible limitation to the values that could take 
the j3 parameter. 

4.2 7= 1/2 and u = : 

G oc A^k~ 2 ~^r 2+ ^ 



c oc AF^kF^r^^m 
A oc AZ^kF^r^ 

Si p = G oc r I, c oc t" 1 / 3 , A oc t" 4 / 3 

c oc t -1 / 3 ,this result also it is obtained by Barrow (see 0) and Petit (0) 
but not by Troiskii (0). The rest of the quantities coincides with the model 
presented by Petit, except for energy density since Petit considers that the 
mass also should vary. 



p oc k 2 t 2 p oc t 



-2 



/ oc A^ k n 3(1 " /3) tw=m if (3 = o => f oc t 2/3 
this model with /3 = is very similar to a FRW with matter predominance. 



5 Conclusions. 

We have solved through D.A. a flat model i.e. the sectional curvature of 
the 3-space is zero, homogeneous and isotropic i.e. we admitted symmetries 
type FRW. The energy-momentum tensor is described by a fluid with bulk 
viscosity in the one which furthermore we envisage mechanisms so much of 
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creation of matter as of entropy and in the one which the classics " constants" 
G, c and A are considered as variable. 

The envisaged cases here show the following behavior for such " constants" 

Goer 1 c oc r 1/2 A oc r 1 

for a model with (7 = 1/2, u — 1/3 y f3 — 0) and 

Goer 2 / 3 ccxr 1 / 3 Acxr 4 / 3 

for a model with (7 = 1/2, u; = y /? = 0) an equal behavior for the rest of 
the quantitites that are observed it in the classic FRW. 

Several problems have emerged during the development of the article. 

1. Even though we have supposed the classic constants of the physics as 
variable have arisen (and result us indispensable) two constants, fc 7 
and Aj,/3 that have clear physical meaning and without those which we 
cannot arrive to solve our model. Will be these two characteristic " 
constants " of our model also universal? 

2. For the calculation of the thermodynamic quantities as the temperature 
and the entropy we have used the "constant" k B that we have supposed 
constant and the "constant" of radiation a that also we have supposed 
constant, but if a is constant (observe that a oc (k%/c 3 h 3 ^) then the 
only one possibility that we have is to make Ti <x. c~ l that is equals to 
say that h oc t 1 ^ 2 or h oc t 1 ^ depending on the model. 

3. If we abandon the characteristic value of 7 = 1/2 we prove without 
difficulty that for 7 > 1/2 the "constant" G varies in a proportional 
way to the time instead of inversely proportional to the time as have 
obtained here. But we would need of some evidence or physical rigorous 
reasoning to take similar values of 7. The only one possibility to obtain 
G and c constant with 7 = 1/2 is imposing a physically unrealistic 
condition u = —1/3. 
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